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^ ■ Abstract 
T— I ■ 

' Data on pp — t- ?7'(958)7r'^7r'^ are presented at nine p momenta from 600 to 1940 MeV/c. 

^ ■ Strong S-wave production of f2{1270)r]' is observed, requiring a J^*^ = 2 ^ resonance with 

Ph! mass M = 2248 ± 20 MeV, T = 280 ± 20 MeV. 

0) , 

The first data are presented on pp — )■ ti'tt^tc^ in flight. These data were taken with the Crystal 
Barrel detector at LEAR. They are part of an extensive study of the J = 0, C = +1 system in 
several channels. Data have been reported earlier on tt^'tt'' [1], 7777 and rjr]' [2], and ?77r°7r° [3]. A 

1) ■ comparison will be made here specifically with the rjn'^iT^ data, and with a combined amplitude 
p • analysis of all the earlier data [4]. 

^ . The experimental set-up has been reported in detail [5]. A p beam from LEAR interacts 

in a liquid hydrogen target 4.4 cm long at the centre of the detector. Incident p are counted 
^ I by a coincidence between a scintillator of 5 mm diameter and a small multiwire proportional 

>■ ■ chamber, both positioned ~ 5 cm upstream of the target. Two veto counters 20 cm downstream 

^ ■ of the target provide a trigger for interactions. The target is surrounded over 98% of the solid 

Q . angle by a multiwire proportional chamber and a silicon vertex detector. These provide an 

^ ! on-line trigger for neutral final states. With a p beam of ~ 2 x 10^/s, the trigger rate is ~ 60/s. 

^ I The present channel is studied in IO7 events, where 77' — )■ rjiT^n^, rj — )■ 77. Photons are 

^ ■ detected with high efficiency down to 20 MeV in a barrel of 1380 Csl crystals covering 98% of 

^ ■ the solid angle; the geometry is such that crystals point towards the target. The crystals have 

>• . a length of 16 radiation lengths and provide an angular resolution of ±20 mrad in azimuth and 

^ ! polar angle. The energy resolution is given by AE/E = 2.5%/-E(GeV)^/^. 

^ I The general procedures for event reconstruction and selection have been described in several 

earlier publications, of which the most detailed concern the study of n^ir^, rjri and 7777' final states 
[1,2]. A Monte Carlo simulation of the detector is used to assess the efficiency for reconstruction 
of the Tc^TT^T]' final state and the levels of background from competing channels. 

Events are first submitted to a kinematic fit to pp IO7, requiring a confidence level > 5%. 
The best kinematic fit to pp — )■ r]4n^ is then selected, again with confidence level > 5%. At 
this step, the main background to rjAn^ comes from Svr^ events. This background is suppressed 
strongly by rejecting any event passing a kinematic fit to bir^ with confidence level > 1% (or 
0.1% at 600 MeV/c, where the background is more severe). Finally, those few events are rejected 
which fit rjTj'iTi^ with confidence level better than rjAn^. 

Fig. 1 illustrates at four beam momenta the t^tttt mass distribution of surviving events in 
the mass range around the 77'. There is a clear rj' signal, agreeing in mass within < 4 MeV 
with the standard value at all momenta. It is superposed on a smooth background, whose 
magnitude is largest at low beam momenta. The Monte Carlo simulation estimates that the 
background comes approximately equally from 3 sources: (i) Svr^ events, (ii) ujAt[^ , [u — ?■ -k^'-j) 
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Figure 1: The distribution of M^{f]7r7r) at four beam momenta indicated by numerical values in 
MeV/c. The shaded areas show selected signal events in the rj' peak and those used for sideband 
subtraction. 



after losing one photon, and (iii) rjAiT^ without an 1]'. The predicted background agrees with that 
observed (within 10% of the prediction). Tighter cuts do not improve the signal/ background 
ratio significantly, but simply cause loss of events. 

Signal events are selected from the peak region of the r]' by adjusting a mass cut around the 
peak at every individual momentum so as to optimise the signal/background ratio. Very rarely, 
two events fall within the window; in this case the one closer to the 1]' is accepted. Statistics 
of the data selection are shown in Table 1. In the maximum likelihood fit used for amplitude 
analysis, sidebins events shown shaded in Fig. 1 are used to subtract the background. The 
areas of sidebins are chosen so that each covers twice the range of mass squared which is used 
to select f]' events; in this way, statistical errors on the background arc small. A technicality is 
that the width of the mass cut is varied according to the accuracy with which the r]' mass is 
reconstructed. This is the reason that sidebands have diffuse edges: the width of the sidebin 
likewise varies with the width of the rj' mass cut. Technically, the way the subtraction is made 
is to include sidebin events into the fit with a weight —0.25 times that of events selected in the 
signal region. Amplitudes are constructed with tensor expressions using the measured mass of 
each 1]'. 

Fig. 2 shows the Dalitz plots at all momenta for events from the signal region. There 
is an obvious contribution due to /2(1270)?7', appearing at momenta > 1200 MeV/c at the 
lower left edge of the plot. Fig. 3 shows the Dalitz plots for sidebin events. The distribution of 
background is not uniform, but peaks in the corners of the Dalitz plots. This peaking accounts for 
corresponding peaks observed in the corners of the Dalitz plots of Fig. 2. When the subtraction 
is made, the surviving signal outside the /2(1270) peak is nearly uniform within the available 
statistics. At 1940 MeV/c, there is also some weak /2(1270) in the background; we have checked 
that this is not due to rj'n^n^ signal spilling into the mass ranges used for the sidebins. 
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Figure 2: Dalitz plots at all beam momenta for events from the signal region of Fig. 1. Numerical 
values indicate beam momenta in Mev/c. 
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Figure 3: Dalitz plots for events from the sidebin regions of Fig. 1. 
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600 


180 


61 


119 


2.90 


900 


1017 


399 


618 


4.61 


1050 


831 


257 


574 


5.76 


1200 


2770 


852 


1918 


6.33 


1350 


2296 


595 


1701 


5.92 


1525 


1416 


381 


1035 


5.06 


1642 


1530 


330 


1200 


4.72 


1800 


1503 


325 


1178 


4.57 


1940 


1063 


240 


823 


4.34 



Table 1: Numbers of selected events, estimated background (BG), true signal, and reconstrucion 
efficiency e as a function of beam momentum. 



There is no indication for the presence of 02(1320) — )■ rji^. The expected contribution may 
be predicted from fits which have been made to a2(1320)7r in rjTT^TT^ data [3]. The predicted 
contribution is only ~ 3% of rj'n'^n^, because of the small (0.53%) branching fraction of 02(1320) 
to Tj'TT. This contribution is included in the amplitude analysis using amplitudes fitted to the 
riim data, but is so small as to have negligible effect on conclusions. Fig. 4 shows projections at 
two beam momenta on to masses of tttt and nr]; the latter is featureless. The histograms show 
results of the maximum likelihood fit described below. 
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Figure 4: Projections on to M(7r7r) and M{r]'7r) at beam momenta of 1050 and 1800 MeV/c; in 
all background subtraction is made using sidebins. Histograms show the fit compared 

with data. 



We now turn to physics results. Data points on Fig. 5(a) show the integrated 77'7r°7r° cross 
section after background subtraction and after scaling to allow for all other unobserved decay 
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modes of r( ^ rj and 7r°. There is a peak around 2230 MeV, which is the nominal threshold for 
/2(1270)r7'(958). The absolute normalisation is obtained using beam counts, target length and 
density, and correcting the observed number of signal events for the reconstruction efficiency 
shown in Table 1. A correction is applied for observed dependence of the cross section on beam 
rate, as described in detail in Ref. [1]. 

The amplitude analysis is made using (a) S and P-waves for at]', where a stands for the tttt 
S-wave amplitude, for which we use the parametrisation of Zou and Bugg [6], (b) S and P-waves 
for /2(1270)?7', and (c) a small, almost negligible contribution from ^Pi — ?> fo{975)ri', which 
helps marginally in fitting the tttt mass distribution at the lowest three beam momenta. It is 
to be expected that higher partial waves for f2r}' will be suppressed strongly by the centrifugal 
barrier in the final states. Contributions from f2i]' D-waves have been tried in the fit, but are 
not required; indeed, the P-wave contribution is quite small. Likewise, arj' contributions with 
L > 2 are negligible. 

We shall present amplitudes for /2(1270)r]' in partial waves (J^*" = 2""'"), and ^F2 
(2++), ^Pi (1++) and ^Fg (3++); they will be compared with /2(1270)77 observed in rjnn data 
[3,4] . These two channels are related by the composition of the rj' and rj in terms of strange and 
non-strange quarks: 

I ^ uu-\- dd 

Iry > ~ 0.8 p 0.6ss, (1) 

v2 

, , uu + dd 

IV > ~ 0.6 ^ + 0.8ss. (2) 

v2 

The coefficients 0.8 and 0.6 are derived from the well known pseudo-scalar mixing angle [7]. Our 
earlier analysis of pp — )■ tt^tt''", 7r°7r°, rji] and rji]' [8] finds that almost all s-channel resonances 
produced in pp interactions are consistent with small mixing angles < 15° between [uu + dd] / \/2 
and ss. The naive prediction is therefore that amplitudes a for pp /2(1270)r;' and pp 
/2(1270)77 will be related by 

a{f2r]') ~ 0.75a(/2?7). (3) 

The peak in the full curve of Fig. 5(a) requires a resonance in f2ri' close to the mass of the 
peak. However, the mass spectrum from a simple resonance will be pushed upwards by the 
rapidly increasing phase space for the final state /2?7'. This effect is visible in the dotted curve 
of Fig. 5(c), which shows the resonance contribution to f2ri' fitted to ?72(2248); this curve peaks 
above 2300 MeV because of the increasing phase space. In order to reproduce the integrated cross 
section of Fig. 5(a), the amplitude analysis requires a strong interfering background peaking 
below threshold. The interference is constructive at low masses, and is required to give a large 
/2?7' cross section there, despite the limited phase space. Above the peak at 2230 MeV, the 
interference becomes destructive, and cuts off the f2r]' cross section on the upper side of the 
resonance. 

The motivation for including this background contribution at low /2?7' masses arises from 
the new combined analysis [4] of TyTTTr data, together with those on pp — >■ tt^tt"*", 7r°7r°, r]r] and 
7777'. Results for r^vrvr from that analysis are shown in Fig. 5(b). That analysis requires a 
resonance at 2267 ± 14 MeV. It appears there most clearly in /2(1270)?7 with L = 2 in the final 
state, shown by the chain curve in Fig.5(b). However, for the dominant f2r] L = channel, 
what one observes is a strong peak near 2 GeV, shown by the full curve. This comes mostly 
from 772(1860), but partly from 772(2030) reported in an analysis of data on pp — > rjn^n^n^ [9]. 
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Mass (GeV) 

Figure 5: (a) Points with errors show the integrated cross sections for the final state jyV^Tr^, 
after correction for backgrounds and for all decay modes of r^', rj and 7r°; the full curve shows 
the fit from the amplitude analysis; the dashed curve shows the rjTnr cross section from Ref. 
[3], multiplied by the SU(3) factor (0.75)^; the dotted curve shows the ar]' contribution; (b) 
the full curve shows the cross section for f2r] fitted to t^tttt data; the dotted curve shows the 
contribution to TyTTTr from 772(2248) alone and the dashed curve that from 772(1860) + ?72(2030); 
the chain curve shows the intensity fitted to 72(1270)77 with L = 2 in Ref. [4]; (c) as (b) for 7277'; 
(d) Argand diagram for the f2i]' S-wave amplitude; crosses mark beam momenta; (e) Intensities 
of contributions to f2i]' from (full curve), ^Pi (chain curve), (dotted) and ^^2 (dashed); 
(f) intensities of contributions to at]' from 0~ (full curve), 1+ (dashed) and l"*" /o (975)77' 
(dotted). 
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The intensities of contributions to the f2r) channel are shown in Fig. 5(b) from (i) all ri2{2248) 
contributions (dotted curve) and (ii) the coherent sum of 772(1860) and 772(2080) (dashed curve); 
the latter two resonances are not well resolved by the rinn data, because they lie close together 
near the pp threshold. The contribution from 772(2248) interferes destructively with 772(1860) 
and 772(2030), so as to cut off the full curve at high masses. 

In present data, the width of the 772(2248) is well determined by the width of the peak in Fig. 
5(a): r = 280 ± 20 MeV. This determination is superior to that in rjTnr data: 290 ± 50 MeV. 
The mass is somewhat less well determined, since the interference with the tails of the lower 
resonances may shift the peak by an amount which is sensitive to their widths. Using the best 
estimates for the widths from Ref . [8] , the mass from the present data is M = 2248 ± 20 MeV, 
in reasonable agreement with the value derived from tjtttt data: 2267 ± 14 MeV. The Argand 
diagram for the 7277' S-wave amplitude is shown in Fig. 5(d). 

A striking feature of the 7277' signal is its large magnitude. The dashed curve on Fig. 5(a) 
shows the complete integrated rjn^n^ cross section, multiplied by (0.75)^ to allow for the expected 
inhibition of 77' with respect to rj. It is surprising that the /277' signal is nearly as strong as the 
dashed curve, bearing in mind the difference in available phase space for 7277' and 7277. The peak 
in the 77V7r cross section (full curve) is much larger than the small peak observed at the same 
mass in the rj-n-R cross section. Likewise, the S-wave peak due to 772(2248) — ?> 7277', shown by the 
dotted curve in Fig. 5(c), is considerably stronger than that in /27/ in Fig. 5(b). If one takes into 
account the available phase space for 7277' and /277, the couphng constant for 772(2248) — > 7277' 
relative to that in /277 is stronger than predicted by equn. (3) by a factor 5.2 in amplitude. 

Vandcrmeulcn has remarked that pp annihilation usually favours high mass final states [10]. 
This may be understood as a form factor effect, arising from the sizes of the participating 
states. In present data, the final state /27/' has very low momentum. However, in the process 
772(2248) — )■ /27^, the momentum q in the final state is ~ 635 MeV/c. The factor 5.2 would 
require a form factor exp— (4.1^^) in amplitude, with q in GeV/c; if this arises from a source 
having a Gaussian distribution in r, the form factor takes the well known form exp —{q'^R'^/6), 
and requires a radius of interaction R = 0.98 fm. Such a form factor is surprisingly strong. For 
comparison, the Vandermculcn form factor approximates to exp — (1.5g^). 

A possibility is that 772(2248) is an ss state. However, strong production from pp is unlikely 
and in disagreement with results for tttt, 7777 and 7777' [4]. 

The strong sub-threshold contribution to the f2f]' S-wave is intriguing. A variety of explana- 
tions are possible, of which we mention one. In Ref. [9], evidence has been presented for three 772 
resonances in a mass range where only two are likely to be qq. Of these, 7/2(1860) is a candidate 
for a hybrid, because of its strong decay to /27/, despite limited phase space. If that conjecture is 
correct, it should be accompanied by an ssg partner at about 2100 MeV. Such an ssg hybrid is 
expected to decay strongly to /2(1525)77' and 72(1270)77'. If it mixes into neighbouring qq states, 
it could help to explain the anomalously strong f2f]' signal observed here. 

We now consider other partial waves. The present data require a small but significant P-wave 
/277' contribution. This could arise from initial pp states ^Pi, ^P2, ^-^2 or '^^3. The amplitude 
analysis of Ref. [4] requires all of these contributions in t/tttt data with a 3+ resonance at 2303 
MeV, a 1+ resonance at 2310 ±60 MeV and 2+ resonances at 2240 and 2293 MeV. A good fit to 
present data may be obtained by fixing the relative magnitudes and phases of these partial waves 
from the fit to rjTT^n^ data. The absolute magnitude of the P-wave contribution is sensitive to 
the radius chosen for the Blatt-Weisskopf centrifugal barrier. This radius is therefore adjusted 
to give the best fit to the data, with the reasonable result 0.8 fm. 
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Figure 6: Production angular distributions for /2?7' at (a) 1525 and (b) 1940 MeV/c for M^tttt) > 
1.1 GeV; also for ar]' at (c) 900 and (d) 1800 MeV/c for M{Tnr) < 1 GeV. Points with errors 
show data, uncorrected for acceptance; histograms show the maximum likelihood fit. 



The magnitudes of the contributions are then 3.5% for ^F^, 3.2% for ^Pi, 3.2% for the 2"*" 
resonance at 2240 MeV and 1.0% for the 2"*" resonance at 2293 MeV; in the latter two, the ratios 
of amphtudes for and are taken from Ref . [4] . Without these amphtudes, log likehhood 
of the fit to ri'n^n^ is worse by 142 for only one parameter fitted to the overall magnitude; so 
the P-wavc contribution is highly significant. [Our definition of log likelihood is such that it 
a change of 0.5 corresponds to one standard deviation change in one variable]. If instead the 
magnitudes and phases of these amplitudes are fitted freely, the fit changes very little. It is not 
possible from the present data to separate and ^F2, which need to be constrained in relative 
magnitude as determined in Ref. [4]. With this constraint, the freely fitted intensities are 3.9% 
for ^Fa, 4.7% for ^Pi and 3.9% for 2+, close to the contrained fit. 

Figs. 6(a) and (b) show angular distributions for production of rj'f 2{1270) in the mass range 
> 1 . 1 GeV in terms of the centre of mass angle 6 of the tj' The distributions are uncorrected for 
acceptance, which is included in the maximum likelihood fit shown by the histograms. At high 
beam momenta, the acceptance for rj' falls in the forward direction, where the separation of its 
decay products becomes less efficient. A check on the reconstruction procedure is that angular 
distributions are symmetric forward-backward in the centre of mass system within errors, after 
correction for acceptance; this symmetry is required by charge conjugation invariance. 

We now turn to the contributions from the broad arj' channel. Prom present data, the only 
firm conclusion which may be drawn is that contributions from both ^5*0 and ^Pi initial states 
are required. At all momenta from 900 MeV/c upwards, the data require angular distributions 
of the form A + Bcos^Orf', as shown in Figs. 6(c) and (d). The t^tttt data have been fitted 
including 0~ and 1"*" resonances. Present data are fitted well by the same resonances. However, 
statistics are not sufficient to provide clear evidence of these resonances in present data. Fig. 6 
shows that the fit to data is adequate. 
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In summary, the main feature of the r/Vvr^ data is a peak at 2230 MeV, requring a dominant 
contribution from the /2(1270)?7' S-wave. The data require a 2~~^ resonance with mass 2248 ±20 
McV and width T = 280 ± 20 MeV; this result is closely consistent with an ri2{2267) resonance 
observed in rjim data. The /2?7' S-wave amplitude is surprisingly strong compared with that for 
/277, even allowing for a form factor in the latter. Contributions from /2(1270)77' P-states are 
consistent with the amplitude analysis of the rjinr data. 
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